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The annihilation of light dark matter was recently advocated as a possible explanation of the large 
positron injection rate at the Galactic center deduced from observations by the SPI spectrometer 
aboard INTEGRAL. The modelling of internal Bremsstrahlung and in-flight annihilation radiations 
associated to this process drastically reduced the mass range of this light dark matter particle. 
We estimate critically the various energy losses and radiations involved in the propagation of the 
positron before its annihilation — in-flight or at rest. 

Using a simple model with mono-energetic positrons injected and confined to the Galactic bulge, 
we compute energy losses and gamma-ray radiations caused by ionization, Bremsstrahlung inter- 
actions as well as in-flight and at rest annihilation and compare these predictions to the available 
observations, for various injection energies. 

Confronting the predictions with observations by the CGRO /EGRET, CGRO/COMPTEL, IN- 
TEGRAL/SPI and INTEGRAL /IBIS /ISGRI instruments, we deduce a mass upper bound of 3 to 
7.5 MeV/c 2 for the hypothetical light dark matter particle. The most restrictive limit is in agreement 
with the value previously found by Beacom and Yiiksel and was obtained under similar assumptions, 
while the 7.5 MeV/c 2 value corresponds to more conservative choices and to a partially ionized prop- 
agation medium. We stress how the limit depends on the degree of ionization of the propagation 
medium and how its precision could be improved by a better appraisal of data uncertainties. 

PACS numbers: 95.30.Cq,95.35.+d,95.55.Ka,98.35.-a,98.38.Am,98.70.-f,98.70.Rz,98.70.Sa 



I. INTRODUCTION 

Observed by various instruments since the 70's the 
511 keV gamma-ray line emission from the Galactic cen- 
ter has now been studied for almost four years with the 
SPI spectrometer @ aboard the INTEGRAL satellite. 
While this emission is definitely the result of electron- 
positron annihilations in the central region of the Galaxy, 
the origin of these positrons remains obscure. 

What is clear, however, from the spectral analysis 
of this line, is that most of the positrons are set to 
rest, form positronium and annihilate in the Galactic 
bulge (hereafter GB) H H Q. Thus, there should 
be material in the bulge to slow down these leptons, 
as well as a significant magnetic field to confine them, 
which is an important information since the interstel- 
lar medium in this region is very difficult to observe. 
As far as positrons are concerned, the main character- 
istic to explain is their very high injection rate, i.e. 
~ 1.4 X 10 43 e + s _1 , as deduced from the 511 keV line flux 
of (1.07 ± 0.03) x 10~ 3 phcm" 2 s" 1 measured by INTE- 
GRAL/SPI .4]. Though previous instruments like SMM, 
TGRS , CGRO/OSSE gave essentially the same require- 
ment |l| , it is surprising that nobody had remarked that 
type la supernova?, a long time considered the main 
source of GB positrons, fall short explaining this high 
injection rate [a, Sill- Fortunately, this drawback is now 
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taken seriously. 

Another clue to the origin of these GB positrons is the 
large extent of the emitting region, 8 ° FWHM Q , in- 
dicating that either they emanate from a single central 
source or a population of sources densely packed in the 
center of the Galaxy and diffuse to fill the whole bulge 
region, or there is a source population with an exten- 
sion similar or greater than the GB. Indeed, the 511 keV 
emission region traces where positrons annihilate, and 
is therefore related to the distribution of the gas in the 
Galactic bulge. 

Potential astrophysical sources of positrons are numer- 
ous (SNIa, hypernovae, gamma-ray bursts, cosmic rays, 
low mass X-ray binaries, millisecond pulsars,. . . ) and 
we will not discuss them in detail here. None of them 
qualifies to account for both the flux and morphology 
of the 511 keV emission. Some candidates are excluded 
on the basis of the disk over bulge ratio [j| and oth- 
ers like millisecond pulsars |10| and gamma-ray bursts 
[Til are shaded off, the first because the electrons 
and positrons would be injected at such energies that 
their Bremsstrahlung emission would be excessive, and 
the second because their rate of appearance is insufficient 
to supply positrons at the required rate 0, Q, 0] . 

In front of this situation, two rather daring hypotheses 
have been put forward : 

- the sources belong to the established list but there 
is something important not taken into account 
in the astrophysical scenario. Positrons from the 
Galactic disk can be transferred to the bulge where 
they annihilate [jj. However, the transport mecha- 
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nism from the disk to the center is poorly known, 
depending critically on the topology and strength of 
the large scale magnetic field and on the diffusion 
coefficient of low energy positrons in the relevant 
magnetic field 0, |n| ; 



- new sources or mechanisms are required, which are 
not in the catalog. We will consider this second 
proposal, and more specifically the ligh t dark mat- 
ter (hereafter LDM) scenario H H H 13 El 
Other exotica from the particle physics community 
include balls [23], relic particl es pll| . decaying 
axinos [22j , primordial black holes [23j , color super- 
conducting dark matter [24|, superconducting cos- 
mic strings pHj , dark energy stars (2|| , and moduli 
decays |22|; they will not be discussed here. 

Surprisingly, though the 511 keV emission from the GB 
was first detected more than thirty years ago, it is only 
with the SPI observations that a great deal of imagina- 
tion and activity has been triggered among the astropar- 
ticle community. 

In this contribution, we gather 7-ray observations per- 
formed by the CGRO (COMPTEL, EGRET) and INTE- 
GRAL (IBIS/ISGRI, SPI) missions and see what con- 
straints they set on the LDM candidate. A note of warn- 
ing however: since we are interested in diffuse radiation, 
the delicate problems of unresolved sources, particularly 
acute at low energies, and of alternate diffuse sources 
should always be kept in mind. 

We do not focus on the theoretical morphology of the 
LDM emission [2£|, which is a thorny subject. It 
first depends on the dark matter distribution, which is 
not well established, varies considerably from one author 
to the other [see e.g. 12^ . [3(j and could furthermore be 
clumpy. It also depends on the complex spatial prop- 
agation of the leptons resulting from LDM annihilation 
and on the distribution of the ISM matter with which 
they interact, also poorly known. Therefore, we use the 
morphological information we dispose of, i. e. that of the 
radiation from positron annihilations occurring at rest, 
as observed by INTEGRAL in the 511 keV line and the 
ortho-positronium continuum. 

We content ourselves to calculate the total secondary 
emission of electrons and positrons in the course of their 
propagation, the only free parameter being the mass of 
the LDM particle. Energy losses are central to the prob- 
lem, since all cross-sections of interest depend on energy, 
except internal Bremsstrahlung associated to the pair 
production. We then combine the resulting global the- 
oretical spectra with our morphological assumptions to 
confront them to observations. The main difficulty is not 
theoretical, but observational: it will be to separate the 
bulge emission in the continuum from that of the disk, 
due to cosmic rays. 



II. OVERVIEW 

When colliding in the central part of the Galaxy, LDM 
particles %, with mass m x < 200 MeV, are expected 
to annihilate into e _ — e + pairs. The total energy im- 
parted to the pair is equal to twice this mass m x , still 
unknown. Internal Bremsstrahlung in the course of the 
pair production [3ll l32l | produces a continuous spectrum 
independent of the density of the surrounding medium. 
Due to the mass limitations, both electrons and positrons 
are intrinsically of modest energy and their energy losses 
(proportional to the interstellar density) are mainly due 
to Coulomb interactions with the surrounding medium. 
In the course of their propagation, they mainly radiate 
through the (external) Bremsstrahlung process. More- 
over, the released positrons are deemed to annihilate with 
ambient electrons, producing a definite 7-ray signature 
(511 keV line plus a low energy continuum). Positron an- 
nihilation can occur in flight |33j — in 26% of the cases at 
most, for m x = 200 MeV — or at rest — mostly through 
positronium formation. Thus, through the combination 
of these various mechanisms, a specific 7-ray spectrum, 
with a discrete and continuous component, is produced 
by positrons between their creation and their annihila- 
tion. 

Our aim is to set an upper limit to the mass of the 
LDM particle self-consistently, by computing this spec- 
trum and confronting it to the observations in the rele- 
vant energy range, i.e. between ~ 10 keV and 200 MeV. 
From the particle physics point of view this new con- 
straint, combined with distinct and i ndep endent ones 
coming from big-bang nucleosynthesis |34(, type la su- 
pernovae, the extra-Galactic hard X-ray and soft 7-ray 
background 0, |H , should help finding a candidate with 
the appropriate mass and annihilation cross-section prop- 
erties. 



III. MODEL 

We develop a simple steady-state model adequate for 
the low energies we are considering, making the reason- 
able assumption that electrons and positrons below 200 
MeV are efficiently confined to the GB, where they an- 
nihilate, due to their small gyroradius — a common hy- 
pothesis in this field. We also consider that the only 
positrons present in the GB are produced in situ, at vari- 
ance with Prantzos 0, who speculates that there could 
be a transfer of positrons from the disk to the bulge. This 
daring possibility is very hypothetical since the magnetic 
field of the Galaxy at large scales and the propagation of 
low energy positrons are poorly known, as admitted by 
Prantzos himself. Our problem involves simple physics, 
repeatedly employed for instance in solar flares physics 
f|3(| and references therein) and Galactic physics [37| . 
and various cross-sections that should be carefully cho- 
sen. 
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A. Steady-state equilibrium 



We consider mono-energetic positrons of total energy 
£7j n j injected in the GB at a rate i?i n j. As they travel 
through the interstellar medium, they loose energy to 
the ambient medium and can undergo in-flight annihila- 
tion. Supposing the positrons are confined in the bulge 
and neglecting their spatial distribution, the differential 
number N(E) of positrons with energy E in steady-state 
obeys to the diffusion-loss equation 



d 

dE 



dt J 



-Q{E) , 



(1) 



where 4f denotes the mean energy loss rate of positrons 
and the source term Q(E) = R in jS(E - E ln j) - Qja{E) 
aggregates the mono-energetic injection and the in-flight 
annihilation, which is a function of energy: 



Q IA (E)=n e a IA (E)v(E)N(E), 



(2) 



with n e the electron density and gia the in-flight anni- 
hilation cross-section. Introducing g(E) = ^z^^n^^ 
and G{E) = g(E')dE' = \ogp(E), where p(E) is the 
probability for a positron injected at E- ln j to survive until 
energy E, the number of positrons in the bulge follows 
from Eq.[f]andEl 



N ( E ) = -Mdt eG[E) (™ c2 < E < E ^ 



(3) 



We can then obtain the total rate of positrons undergoing 
in-flight annihilation 



Ria = 



Q IA (E)dE = R inj [l-e G< - mc >}e + s 



(4) 



and deduce the rate of annihilation at rest, assuming it 
equals the rate of positrons reaching complete rest: 



R A = i? inj e G(mc2 ) e+s- 1 . 



B. Energy losses 



(5) 



The total energy lost by positrons during thermaliza- 
tion goes into various interactions - excitation or ioniza- 
tion, Bremsstrahlung, inverse Compton and synchrotron 
effects. As their importance varies with the positron's en- 
ergy, it is mandatory to use refined energy loss rates. For 
the range of injection energies we are interested in (up to 
200 MeV at most), we can limit ourselves to ionization 
and Bremsstrahlung. 

We stress that the ionization degree Xi — nH II jnu 
in the region of propagation, where tih is the hydro- 
gen number density, is of prim e importance for the es- 
timate of Coulomb losses j3j|, since at a given energy 
they are several times weaker in a neutral medium than 
in a plasma. 




FIG. 1: Energy losses of positrons in a 51% ionized medium 
(uh = 0.1 cm -3 ): energy loss rate per unit path length ^ as 
a function of the positron's Lorentz factor 7. Losses through 
Coulomb collisions within an ionized medium (dotted) prevail 
on ionization of the neutral phase (dashed), but both greatly 
exceed Bremsstrahlung losses (plain). 



Thus, the total loss rate sums up three terms: 



dE 
~dt 



dE 1 
~dt 



Off) 



dE 



I.HI 



dt 



-{n H ,Xi 



dE 1 '" 11 
dt 



(n H ,Xi 



(6) 

where the energy lost by a positron due to ionization of 
neutral matter is taken from j3^ while the loss rate for 
excitation of the ionized component follows |4jj . The en- 
ergy loss caused by Bremsstrahlung interactions is given 
by Eq. 4BN in [41| . The relative amplitude of each com- 
ponent is illustrated in Fig.^ which shows energy losses 
4^ per unit path length. 



C. Radiation spectra 



1. Positron radiation 

The total radiation spectrum due to positrons is the 
sum of five components: the 511 keV line from e + -e _ 
annihilation at rest and from para-positronium decay 
into two gamma rays (i), the ortho-positronium spec- 
trum from its 37 decay (ii), the continuous spectrum 
from e + -e~ in-flight annihilation into 27 (iii), the ra- 
diation spectra associated with Bremsstrahlung energy 
losses (iv) and with the internal Bremsstrahlung taking 
place in the annihilation process of the dark matter par- 
ticle (v). 

The flux from 27 annihilation at rest F 21 is tuned 
via the injection rate i?i n j to the value of (1.07 ± 
0.03) 10" 3 ph/cm 2 /s deduced from INTEGRAL obser- 
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FIG. 2: Differential in-flight annihilation cross-section -given 
for instance by Aharonian and Atoyan |43l Eq. 5, mis- 
print corrected]- as a function of the reduced photon energy, 
for various positron energies. Still energetic positrons annihi- 
late with thermal electrons into two photons that share the 
total available energy, with a minimum energy of m e c 2 /2 by 
photon. 



vat ions Q 



_p2 7 



(l-fps) + f Ps /i 



4nd 2 GG 



Ra, 



(7) 



using the associated positronium fraction value 
f Ps = 0.967 ± 0.022 0] and the relationship © be- 
tween the initial positron injection rate and the rate Ra 
of positron annihilation at rest. 

The 37 ortho-positronium decay spectrum follows 



dE-, 



' 4nd 2 GC 



(8) 



where the spectral distribution ^ L was derived by Ore 
& Powell H2 in ll949i 

Convolving the differential in-flight annihilation cross- 
section (Fig. EJ with the positron distribution yields the 
in-flight annihilation component: 



dF 



I A 



1 



da 



1 A 



4ird GC JEmin(E^) 



-{E 1 ,E)v{E)N{E)dE . 

The external Bremsstrahlung s pect rum is derived in a 
similar way using Koch & Motz g^'s cross-section 3BN. 
Dependence on the electron number density n e cancels 
out as both the differential emissivity n e -j^-v and the 
energy losses in the positron spectrum N(E) (see Eq. O 
are directly proportional to n e . Finally, the internal 
Bremsstrahlung spectrum depends only on the positron 



dF 



IB 



1 da 



IB 



dEry And GC a to t dE 1 



(10) 



where we use the differential cross-section first derived 
by Beacom et al. jU- Another cross-section has been 
derived by Boehm & U wer since then, which appears 
to be larger. 



2. Electron radiation. 

In addition to the five positron radiation components 
comes the Bremsstrahlung spectrum from the mono- 
energetic electrons injected together with positrons dur- 
ing the dark matter annihilation process, whose equilib- 
rium spectral distribution follows 



N e -(E) = - 



Ru 



dE/dt 



(mc 2 < E < E h 



(11) 



D. Discussion 



Fig. |3| displays the various components of the radiation 
emitted by positrons and electrons from dark matter an- 
nihilation for a specific mono-energetic injection. While 
external Bremsstrahlung becomes the main component 
at very low energies, it is negligible compared to both 
in flight annihilation and internal Bremsstrahlung above 
1 MeV, which we will show to be the most constraining 
energy domain. Fig. 0] shows the evolution of the total 
spectrum with increasing injection energies and the in- 
fluence of the ionization fraction. 

Using a multi-component model of the interstellar 
medi um (molecu lar, cold, warm neutral, warm ionized, 
hot). I Jean et alJ Q showed, through their analysis of the 
511 keV line and the ortho-positronium spectrum, that 
51^2% of the rest annihilation takes place in the warm 
ionized phase of the ISM. We can assume that the mean 
degree of ionization of the medium where the thermaliz- 
ing positrons propagate ranges between and 51%. We 
consider both extreme cases; the latter one enables us to 
release constraints on the mass of the dark matter parti- 
cle. 



IV. CONFRONTATION TO OBSERVATIONS 



In section ITlIl we computed the total diffuse 7-ray flux 
caused by the production, interactions and annihilation 
of positrons and electrons coming from the annihilation of 
a LDM particle (hereafter the LDM flux). We now would 
like to use observations by INTEGRAL (ISGPJ and SPI) 
and CGRO (COMPTEL and EGRET) to derive some 
constraints on the mass of this particle. The computed 
flux stands for the GB, in which we assume that positrons 
(and associated electrons) are produced and trapped. 
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FIG. 3: Light dark matter (_Ei n j=10 MeV) radiation spectra: 
positron annihilation spectra (in-flight: IA; 27 at rest; 37 at 
rest), positron energy loss spectra (internal Bremsstrahlung: 
IB; Bremsstrahlung: B e + ), electron Bremsstrahlung spec- 
trum (B e~) and total spectrum (plain). 
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FIG. 4: Total dark matter radiation spectrum for increas- 
ing injection energies Ei n j (1, 3, 5, 10, 50 and 100 MeV). 
In a partially ionized medium (red dotted curves), higher 
energy losses through ionization lead to lower external 
Bremsstrahlung fluxes than in a completely neutral medium 
(black solid curves). 



Therefore, we would like to show, for each possible 
injection energy of the positrons : 

(i) whether the available data should and do display 
an excess in the bulge region with respect to the 
surrounding Galactic plane; 

(ii) whether the total predicted flux exceeds the ob- 
served one, taking other gamma-ray diffuse com- 
ponents into account. The latter include extra- 



Galactic diffuse emission, diffuse emission from the 
interaction of Galactic cosmic rays (CRs) with the 
interstellar medium and a component due to - 
still — unresolved point sources. 



A. Observational data 

We consider the available measurements for diffuse 
emission in the 20 keV to 200 MeV energy range (data 
points in Fig. |S] and 0) . 

Analyzing the INTEGRAL /IBIS /ISGRI data, 
Krivonos et al. |44[ computed the Galactic ridge X-ray 
spectrum in the 17 to 200 keV range using the imager 
IBIS as a collimated telescope. They derived the residual 
emission from the total field of view remaining after sub- 
traction of the instrumental background and of resolved 
point sources. Considering the ~ 19 width of the 50% 
coded field of view, we will be able to compare the fluxes 
they derived to models of a 20 ° wide region around 
the Galactic center. Below ~ 60 keV, there remains 
a residual flux unexplained by resolved point sources. 
Whether this flux is due to truly diffuse emission or 
can be attributed to a population of unresolved sources, 
like Galactic pulsars, is still matter of debate. Above 
60 keV, no residual flux is detected but 2a upper limits 
on the possible unresolved emission have been set for 
three energy bands up to ~ 190 keV. 

For the diffuse emission observed by the INTE- 
GRAL/SPI spectrometer, we use the spectrum obtained 
by Strong et al. [1| in the 18 keV to 1.018 MeV en- 
ergy rang e. Note that it has also been studied indepen- 
dently byjBouche^OD Oil w ho derived similar conclu- 
sions. IStrong et alJ proceeded by fitting the instrument's 
data with a multi-component imaging model including 
point sources detected by the spectrometer and a set of 
ten maps with various diffuse morphologies. The sum 
of the ten diffuse components represents the map of the 
emission seen by SPI as diffuse. Strong et al. published 
([45l|. Tab. B.4) the fluxes of this map integrated over a 
20 ° x20 ° region at the Galactic center, shown in Fig.JS] 
and [7| The associated uncertainties do not include sys- 
tematic effects. We will compare these fluxes to a model 
of emission in a region of the same size, as a more de- 
tailed morphology of the SPI emission inside this region 
is not available yet apart from the morphology studies 
of the para-positronium and 511 keV annihilation com- 
ponents. We notice a discrepancy between the SPI and 
ISGRI diffuse spectra. This gap could be explained by 
the imager's better ability to detect faint point sources, 
considered as diffuse by SPI. However, there might be a 
need for a more profound analysis. 

For CGRO/COMPTEL, we will start from the 
intensity sky maps of diffuse emission produced by 
IStrong et alt^ in the 1-3, 3-10 and 10-30 MeV energy 
bands. These maps are available with a resolution of 1 ° . 
After subtracting a zero level deduced from high lati- 
tudes, we compute total fluxes in regions of variable sizes. 
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Data points in Fig. [S] and correspond to a 20 ° x 20 ° 
region, for mere compatibility with the ISGRI and SPI 
observations. Uncertainties on these maps are largely 
dominated by systematic effects; an estimation based on 
local fluctuations at high latitudes leads to uncertainties 
of at least 30% 0,0] ■ Hence, for all sizes of our integra- 
tion region, we consider the COMPTEL flux uncertainty 
to be 30%. 

We also consider CGRO/EGRET 7-ray measurements 
above 30 MeV ; the fluxes of data points in Fig. Inland 
Hare obtained by subtracting the isotropic extra-Galactic 
component to longitude profiles and integrating them on 
our region of study; associated systematic uncertainties 
amount to ~T5 to 20%. 



B. LDM morphology hypothesis 

Our model does not include spatial diffusion but sim- 
ply assumes that positrons injected in the Galactic bulge 
remain confined to it; therefore, we have no detailed in- 
formation on the morphology of the LDM emission. How- 
ever, analysis of the SPI data showed that the 511 keV 
component is compatible with a 8 ° FWHM Gaussian 
spatial distribution 3] and that the ortho-positronium 
continuum component had a similar spatial distribution 
[5fj| . Thus, we can safely assume that the rest of the LDM 
emission, e.g. radiation emitted during the injection and 
thcrmalization phases, comes from the same region or a 
smaller one. 

When considering a 20 " x20 ° region, we do not need 
to make a strong assumption on the morphology of in- 
flight annihilation and we can use the total LDM flux 
computed in section ITTT1 Conversely, when considering a 
smaller region, we assume an LDM morphology follow- 
ing a Gaussian distribution with an 8 ° FWHM and we 
scale the total LDM flux by the appropriate fraction cor- 
responding to the solid angle subtended by this region. 

In a first step, we use the 20 ° x20 ° region which is the 
smallest one for which data for all four instruments are 
available. In a second step, we only use the COMPTEL 
data, which bring the most constraining information, and 
we consider circular regions of 5, 8 and 14.6 ° of diame- 
ter corresponding to regions covering respectively 24, 50 
and 90% of the LDM flux. Such smaller regions include 
less diffuse background from cosmic-ray interaction and 
enable to enhance the signal to background ratio, leading 
to more stringent constraints. 

C. Full model approach 

The main problem when comparing LDM theoretical 
emission to the data is to estimate the other diffuse com- 
ponents, mainly the emission due to interactions of cos- 
mic rays with the interstellar medium. In this paragraph, 
we adopt the approach consisting in reproducing the ob- 
servations using a complete model. 
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FIG. 5: Comparison between the total diffuse spectrum ex- 
pected from a 20 ° wide region at the Galactic center and the 
diffuse fluxes measured by IBIS (Krivonos et al. 2006 [i^l : 
magenta), SPI (Strong et al. 2005 0; blue), COMPTEL 
(Stronget al. 1999 gj; green) and EGRET (Hunter et al. 
1997 |43 |: cyan). The model curves are the sum of the CR 
interactions component computed bv lStrong et all [4^ (lower 
black line), and the dark matter related component, com- 
puted in § Hill for injection energies i?j n j = 1, 3, 5, 10, 50 
and 100 MeV, in a neutral (black solid curves) or 51% ionized 
medium (red dashed curves). 



We consider the GALPR0P code developed by Strong 
et al. 0j E3 to model the cosmic-ray propagation and 
the associated 7-ray production. The total CR gamma- 
ray emission is composed of three main components : 
Bremsstrahlung scattering, the inverse Compton effect, 
and pion decays. The relative importance of these com- 
ponents varies with energy. 

Fig-Elshows, for a 20 ° x20 ° region, the spectra of the 
total diffuse emission expected. They sum up the radi- 
ation from CRs, obtained by integration of the GALPR0P 
(version 600203a) skymaps, and the total LDM 7-ray 
flux, for various positron injection energies reflecting the 
LDM particle mass. For each injection energy, both 
the completely neutral propagation medium case and the 
51% ionized medium case are illustrated. 

The GALPROP model reliably reproduces the measured 
diffuse spectrum over a wide energy range, from tens of 
MeV to hundreds of GeV. In particular, it reproduces the 
EGRET measurements satisfactorily, mostly through the 
pion decay bump. On the contrary, at lower energies, this 
CR model fails to account by itself for all of the detected 
diffuse emission. Assuming that this model still correctly 
predicts the diffuse emission due to cosmic-ray interac- 
tions below 30 MeV, this leaves room for diffuse emission 
from a different origin or for unresolved point sources. In 
this section, we tentatively adopt the hypothesis that all 
of the discrepancy in the bulge between the data and the 
GALPROP model is due to the sole LDM emission. 

Then, Fig. [5] permits us to derive for each instrument 
constraints on the maximum mass of the LDM particle. 
With a 95% confidence level, in the more stringent case 
of a neutral medium, we exclude LDM particles with 
masses greater than ^70 MeV/c 2 using the ISGRI up- 
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per limits, 200 MeV/c 2 using the SPI flux above the 
511 keV line, 25 MeV/c 2 using the 10-30 MeV COMP- 
TEL flux. For EGRET, the cosmic-ray interaction mod- 
els suffice to explain the measurements; no additional 
diffuse component such as in-flight annihilation or inter- 
nal Bremsstrahlung is required, which sets the maximum 
LDM particle mass to the lower energy boundary of the 
instrument of 30 MeV. 

But the hypothesis adopted that all of the discrepancy 
between the data and the GALPR0P model in the Galactic 
bulge is due to the sole LDM emission suffers from the 
fact that the same discrepancy can also be noticed in the 
Galactic disk, as depicted in Fig.|H| The longitude profile 
of the COMPTEL flux in the 3-10 MeV band is rather 
flat, shows no excess in the bulge region that could be 
identified as a signature of LDM radiation. In order to 
fill the gap between the data and the GALPR0P model in 
the bulge with LDM radiation, we would need to simulta- 
neously explain the outer disk deficit, by adding another 
model component from a totally different origin which 
would be present only outside the Galactic bulge region. 
Such an hypothesis seems unrealistic. 

Additional tuning of the GALPR0P model permits to 
enhance the predicted CR flux in the COMPTEL en- 
ergy range, notably by using a less conventional model 
with a rapid upturn in the CR electron spectrum |52ll53| . 
However, this idea causes energetic problems and has no 
natural motivation : while the origin of the discrepancy 
with the data is still debated, a population of unresolved 
sources is more likely Moreover, such a model would 
explain the COMPTEL observations over the whole lon- 
gitude range. Thus, we should derive our constraints 
simply through the absence of any excess in the bulge 
with respect to the disk in the COMPTEL data. 



D. Empirical approach 

In this paragraph, we drop the idea of using a full 
model approach which, as discussed previously, has the 
drawback of leaving completely out of consideration the 
fact that, outside the energy range of the positronium |50| 
and the 511 keV line no continuum excess has been 
detected in the bulge with respect to the surrounding 
Galactic plane. 

Therefore, taking note of the more disc-like than bulge- 
like morphology suggested by COMPTEL skymaps and 
profiles above 1 MeV, we assume that most of the emis- 
sion detected by both COMPTEL and SPI have their 
origin in either unresolved faint point sources or diffuse 
mechanisms other than those related to LDM particles. 
Assuming a relative uniformity of this 7-ray background 
along the Galactic plane, we derive the mean COMPTEL 
intensity in a wide region (|Z| < 30°, \b\ < 5°) and use it 
as an estimate of background for smaller regions. An ad- 
ditional signal caused by LDM particles is then searched 
as an excess over this background which should remain 
within data uncertainties. This approach, also used by 



iBeacom and Yuksdi p. is illustrated in the right panel 
of Fig. M 

Fig. [7| displays total diffuse radiation spectra. The 
LDM radiation spectra, for various positron injection en- 
ergies, have been added to a powerlaw spectrum based on 
the data itself, at variance with Fig. [5] where the GALPROP 
cosmic-ray induced spectrum was used as a baseline. We 
took the powerlaw derived for the complete SPI unre- 
solved spectrum by Strong et al. |45j, 



^ = (1.39±0.27)10- 3 (-f-)- 1 - 82±a08 cm- 2 s- 1 McV- 1 , 
dE v ; v MeV ; 

(12) 

which is remarkably compatible with the COMPTEL 
data points. This figure still covers a 20 ° x20 ° region to 
display the ISGRI and SPI data points. However, in or- 
der to derive constraints, different sizes of the integration 
region will be taken into account. 

For this purpose, we consider the 1-3, 3-10 and 10- 
30 MeV COMPTEL energy bands, which are the most 
constraining. Fig. [S] illustrates the constraints on the 
LDM particle mass obtained when using an 8 ° diameter 
region. It displays, as a function of the positron injec- 
tion energy and for each energy band, the sum of the total 
COMPTEL flux and of the LDM radiation flux, scaled 
by a factor 50% corresponding to the fraction of the to- 
tal LDM radiation coming from this region under our 
assumption that LDM radiation follows a spatial distri- 
bution identical to that of the observed 511 keV line, i.e. 
an 8 ° FWHM Gaussian. The la and 2a upper values 
on the COMPTEL fluxes are also shown. By excluding 
injection energies yielding a total COMPTEL + LDM 
flux exceeding the 2a COMPTEL upper values, we de- 
rive the maximum mass of the LDM particle at a 95% 
confidence level. The results are given in Tab.QJ together 
with constraints obtained by considering smaller or larger 
regions, for both the case of a completely neutral propa- 
gation medium and a medium with an ionization fraction 
of 51%. 

In the case of a neutral medium, we set a maximum 
of ~ 4 MeV on the positron injection energy with the 
integration region covering the FWHM of the emission. 
When taking a 5 ° diameter, corresponding to 25% of the 
LDM flux, we maximize the signal to background ratio 
and obtain results similar to those ofH eacom and Yiiksell 
|33| : a maximum mass of ~ 3 MeV/c 2 for the light dark 
matter particle. The precise limit depends on the not 
so well known uncertainties on the COMPTEL measure- 
ments. 

The value of the fraction of ionized phase in the 
medium the positrons propagate through has a small 
but not negligible influence on our LDM particle mass 
constraints. By assuming an ionization fraction of 51%, 
those constraints can be released by a factor of approxi- 
mately two (details in Tab.QJ. 
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FIG. 6: Longitude profile of the 3-10 MeV diffuse flux measured by COMPTEL in a latitude band of ±5 ° (histogram). 
Left : the data is compared to a global model — thick line — including the GALPROP model of diffuse emission from cosmic-ray 
interactions, composed of a dominating inverse Compton component — dotted line — and a smaller Bremsstrahlung component 
— dashed line — , to which we added the model of LDM emission — thin solid line — , computed here for an injection energy Ei n , 
of 10 MeV and assumed to have a Gaussian morphology with 8 ° FWHM . As the discrepancy between the COMPTEL data 
and the GALPROP model concerns high longitudes as well as the Galactic center, trying to explain all of the Galactic bulge 
gap through the LDM emission is not reasonable because the outer-bulge disparity would remain unexplained. 
Right : thus, in jTTTjl we dismiss the idea of modeling the observations completely; we acknowledge the relatively flat profile 
of the COMPTEL data in the inner disk region and use the mean intensity as a baseline, requiring that an eventual excess 
caused in the Galactic bulge by LDM radiation remains within data uncertainties. 



Size 


a* [%} 


Energy band 


[*] 


1-3 MeV 


3-10 MeV 


10-30 MeV 


5 





3 


5.5 


12 


51 


7 


6.5 


13 


8 
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12 




51 


9 


7.5 


13.5 


14.6 





10 


8 


14.5 


51 


23 


12 


17 


20x20 





28 


16 


18.5 


51 


60 


27 


23 



TABLE I: Constraints from COMPTEL data on the maxi- 
mum positron injection energy E[ n j [MeV]. The table gives 
the 2ct (95% confidence) maximum energies deduced from the 
1-3, 3-10 and 10-30 MeV COMPTEL measurements respec- 
tively, for four different integration regions (5, 8 and 14.6 ° 
diameter, 20° width), including 24%, 50%, 90% or 99% of 
the LDM emission respectively — under the assumtion that 
positron in-flight annihilation has a morphology similar to 
that of annihilation at rest. 



V. CONCLUSIONS 

Combining all the radiation components related to 
LDM particles — mostly in-flight annihilation and inter- 
nal Bremsstrahlung — we were able to set an upper bound 
of 3 to 7.5 MeV on the mass of these particles, based 
mostly on COMPTEL data. 

A previous study bv lBeacom and Yiiksell p| had con- 
cluded on a 3 MeV/c 2 upper mass. All in all, this work 
shows essential agreement with their result; no incon- 
sistencies were found when making similar assumptions. 
More conservative morphology choices lead to slightly 
larger values but the effect of the ionization fraction of 
the ISM is more noticeable. Arguably, there seems to be a 
tension between the lower limit established through the 
constraint of supernovae explosions mj and the upper 
limit based on in flight annihilation (|33 and this work), 
but a definite exclusion of the LDM hypothesis is pre- 
mature, since the discussion on the uncertainties on the 
measurements in the COMPTEL energy band, which are 
of crucial importance, is rather qualitative. To advance 
in this context, a refined reassessment of both statisti- 
cal and systematic uncertainties on fluxes in COMPTEL 
maps would certainly prove to be useful. Before defi- 
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E [MeV] 

FIG. 7: Comparison between the total diffuse spectrum expected from a 20 ° wide region at the Galactic center and the diffuse 
fluxes measured by IBIS (magenta), SPI (blue) and COMPTEL (green). Instead of taking the predictions from cosmic-ray 
interaction codes as a base level for our LDM model like in Fig. we now use a powerlaw (blue straight line) 7-ray background 
— already displayed in Fig. |S] — compatible with both SPI and COMPTEL data. 



nitely excluding a light dark matter particle of mass less 
than about 10 MeV/c 2 , further observations of the Galac- 
tic bulge region with INTEGRAL/SPI are mandatory to 
get additional morphological information and refined un- 
certainties on diffuse emission in the 511 keV to 1 MeV 
energy domain. 

We stressed the influence of the ionization degree of 
the interstellar medium on the amplitude of the in-flight 
radiation spectrum : taking into account the ionization 
fraction of the propagation region could allow values of 
the mass of the LDM particle larger by up to a factor two. 
However, additional modeling work is required to ascer- 
tain in which phases positrons propagate during their 
thermalization phase. 

The present conclusions were derived under the de- 
manding assumption that light dark matter annihilation 
alone accounts for all of the 511 keV line. A multi-source 
scenario would enable us to reduce the amplitude of the 
LDM 7-ray continuum and to release once again the limit 
on the LDM particle mass. 

This work focuses on the LDM hypothesis; however, 
by leaving the internal Bremsstrahlung component aside, 
the model presented can be applied to the more general 
case of 7-ray radiation from any mono-energetic source 
of positrons. 

Prospects In order to further restrict the LDM parti- 
cle mass range with the INTEGRAL satellite, a refined 
appraisal of the contribution of faint point sources in the 
50-511 keV energy range is required. In this perspective, 
the new census of sources found by the imager above 



100 keV [54[ will be paramount to fine-tune ISGRI up- 
per limits on diffuse emission. Because the spectrome- 
ter does not have the capacity to resolve as many faint 
sources, an increased recognition of ISGRI detections in 
the analysis of SPI data might prove necessary. 

Beyond flux considerations, whose precision is hin- 
dered by our uncertain knowledge of the radiation from 
cosmic-ray interactions and unresolved sources, morphol- 
ogy studies will be a key point to find out whether the 
511 keV to 8 MeV energy range presents a faint excess 
in the GB. 

In order to take all aspects of LDM into account, the 
present model could be extended to include spatial dif- 
fusion of positrons and electrons during thermalization 
and model the differences in the morphologies of the 
LDM emission at various energies. Furthermore, study- 
ing the influence of the in-flight annihilation and inter- 
nal Bremsstrahlung components on the positronium frac- 
tion derived from SPI data would be interesting. 
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10-30 MeV | 

1 10 10 2 

E inj [MeV] 

FIG. 8: Total integrated flux in the 1-3, 3-10 and 10-30 MeV 
COMPTEL energy bands for an 8 ° diameter region around 
the Galactic center. Horizontal lines show the flux — solid 
lines — derived from the mean COMPTEL intensity in a wider 
fraction of the disk (\l\ < 30°, |6| < 5°) and its la — dashed — 
and 2cr — dotted — upper values. The sum of the COMPTEL 
flux and of the LDM flux is shown as a function of the positron 
injection energy -Einj, for a completely neutral propagation 
medium — circled solid lines — and for 51% ionized medium 
— crossed dot-dashed lines. It should be noticed that combin- 
ing the three COMPTEL energy bands instead of considering 
them independently would bring a sensitivity improvement 
and lead to more severe exclusions; for instance, injection at 
Einj=W MeV is excluded by the 1-3 and 3-10 MeV energy 
bands at 3.6a and 5.6a respectively (xi = 0), while the com- 
bination of both bands leads to a 6.9a exclusion. 
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